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ABSTRACT 

K 4-47 is an unusual planetary nebula composed of a compact high-ionization core 
and a pair of low-ionization knots. Long-slit medium-resolution spectra of the knots 
and core are analyzed in this paper. Assuming photoionization from the central star, 
we have derived physical parameters for all the nebular components, and the (icf) 
chemical abundances of the core, which appear similar to Type-I PNe for He and N/0 
but significantly deficient in oxygen. The nebula has been further modelled using both 
photoionization (CLOUDY) and shock (MAPPINGS) codes. From the photoionization 
modelling of the core, we find that both the strong auroral [O ill] 4363A and [N 
li] 5755A emission lines observed and the optical size of the core cannot be accounted 
for if a homogeneus density is adopted. We suggest that a strong density stratification, 
matching the high-density core detected at radio wavelengths and the much lower 
density of the optical core, might solve the problem. From the bow-shock modelling 
of the knots, on the other hand, we find that knots' chemistry is also represented by 
Type-I PN abundances, and that they would move with velocities of 250 - 300kms~^. 
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1 INTRODUCTION 

Planetary nebulae (PNe) are known to possess a variety of 
small-scale structures that are usually in a lower ionization 
state than the main body of the nebulae. The morphological 
and kinematic properties of these low-ionization structures 
(LISs, Gongalves, Corradi & Mampaso, 2001) vary from type 
to type, in the sense that LISs can appear in the form of pairs 
of knots, filaments, jets, and isolated features moving with 
velocities that either do not differ substantially from that of 
the ambient nebula, or instead move supersonically through 
the environment. 

The total number of PNe that are known to possess 
LISs is 55, i.e. about 10% of all the 527 Galactic PNe imaged 
in filters of high- and low-ionization emission lines (Balick 
1987; Schwarz, Corradi & Melnick 1992; Manchado et al. 
1996). The different types of LISs may be easily seen in 
Corradi et al. (1996). 
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At present, the origin of jets and pairs of knots in PNe 
is not completely clear. From the theoretical point of view, 
the principal physical process behind the formation of col- 
limated LISs is the interplay between the stellar AGB and 
post-AGB winds (for single stars) or between stellar and disk 
winds (if the central star is binary). According to the vari- 
ous studies dedicated to this issue (e.g. Garci'a-Segura 1997; 
Garcia-Segura & Lopez 2000; Steffen, Lopez & Lim 2001; 
Blackman, Frank & Welch 2001) jets and knots originated 
by this interplay are predicted to be supersonic, highly col- 
limated and two-sided. In the case of single stars these are 
expected to be produced at the same time as the main PN 
shell, but in the case of a binary star origin they may be 
younger than the host PN. However, important properties 
of these LISs such as density contrasts, the peculiar nitrogen 
abundance and main excitation mechanisms, appear hard to 
explain (see Dwarkadas & Balick 1998; Balick & Frank 2002; 
Gongalves et al. 2003). 

K 4-47 (PN G149.0-I-04.4) is a compact PN that con- 
tains LISs. It is composed by a small, high ionization neb- 
ular core and a pair of low-ionization, high-velocity knots, 
connected to the core by a much fainter low-ionization lane 
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(Corradi et al. 2000; see Figure 1). Gongalves et al. (2001) 
have proposed that the low-ionization lanes and knots of 
K 4-47 are genuine jets. Their morphological and kinema- 
tical properties are explainable if the jets and knots were 
formed by accretion disks, attaining velocities of several hun- 
dred kilometers per second (the main properties of these 
models are also summarized in Gongalves et al. 2001 and 
Balick & Frank 2002). These highly supersonic velocities 
imply that the resulting LISs are likely to be shock-excited. 

K 4-47 is a poorly studied PN, for which statistical 
methods provide an unrealistic wide range of distances, 
for instance, 8.5 kpc (Cahn, Kaler & Stanghellini 1992) or 
26 kpc (van de Steene & Zijlstra 1994). Corradi et al. (2000) 
computed a distance between 3 kpc and 7 kpc assuming that 
the object participates to the ordered rotation of the disk of 
the Galaxy, but note that the relatively large height of K 4- 
47 on the Galactic plane (0.54 kpc for a distance of 7 kpc) 
adds some further uncertainty to this determination. Tajitsu 
& Tamura (1998) estimated a distance of 5.9 kpc using the 
integrated IRAS fluxes under the (crude) assumption of con- 
stant dust mass for all PNe. Lacking of anything better, we 
shall adopt in the following this distance of 5.9 kpc. Lumsden 
et al. (2001) mapped the H2 emission from K 4-47 finding 
that it is excited by shocks. The object also appears in the 
6 cm VLA radio survey of Aaquist & Kwok (1990) showing 
a very compact radio core, with a diameter of 0.25 arcsec, 
and one of the largest brightness temperatures (Ti,=8700 K) 
found in PNe. So far, the properties (luminosity and tem- 
perature) of its central star as well as its nebular (physical 
and chemical) properties are not known. 

In this work, we address the debated issue of the nature 
and origin of high- velocity LISs, through the determination 
of the physical parameters, excitation and chemistry of K 4- 
47. Spectra of the core and the pair of knots are analyzed 
using two different models, which consider the gas to be 
either fully photoionized by the PN central star, or fully 
ionized by shocks. We will show that K 4-47 is particularly 
interesting for this study because, in contrast to most of the 
LISs studied up to now (Dopita 1997; Dwarkadas & Balick 
1998; Miranda et al. 2000; Gongalves 2003; Gongalves et 
al. 2003), its pair of LISs is mainly shock excited. 



2 DATA ACQUISITION AND REDUCTION 

Spectra of K 4-47 were obtained on 2001 August 28 at the 
2.5m Isaac Newton telescope (INT) at the Observatorio del 
Roque de los Muchachos on La Palma (Spain), using the 
Intermediate Dispersion Spectrograph (IDS). The IDS was 
used with the 235 mm camera and the R300V grating, pro- 
viding a spectral coverage from 3650 A to 7000 A with a 
spectral reciprocal dispersion of 3.3 A pix~^ and a resolution 
of 6.9 A. The spatial scale of the instrument was 0".70 pix^^ 
with the TEK5 CCD. Seeing varied from I'M to 1".2. The 
slit width and length were 1.5" and 4', respectively. The slit 
was positioned through the centre of the nebula at P.A. = 
41°, passing through the knots, and the exposure times were 
3x1800 s. 

Bias frames, twilight and tungsten flat-field exposures, 
arcs and exposures of standard stars BD-l-332642, Cyg 
OB2#9, HD19445, and BD-l-254655, were obtained. Spectra 
were reduced and flux calibrated using the standard IRAF 
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Figure 1. [N 11] 6583A image of K 4-47, adapted from Corradi 
et al. (2000). Labels of the selected structures are indicated. 
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Figure 2. Spatial intensity profile of [O l] 6300A, [O 
11] (3726A-I-3729A) and [O in] 5007A emission lines along 
P.A.=41°. 

package for long-slit spectra. Line fluxes were measured se- 
parately for the three nebular regions indicated in Figure 
namely Knotl, Core and Knot2. The observed line fluxes are 
given in Table 1. Errors in the fluxes were calculated taking 
into account the errors in the measurement of the fluxes, as 
well as systematic errors of the flux calibrations, background 
determination, and sky subtraction. The bottom lines of Ta- 
ble 1 give the estimated accuracy of the measured fluxes 
for a range of line fluxes (relative to H/3) in each of the re- 
gions. Absolute H/3 fluxes integrated along the slit in each re- 
gion are as follows: F(H/3)Knoti=2.41; F(H/3)corc=2.76; and 
F(H/3)Kiiot2 = 1.64 (in units of lO"^'^ erg cm"^ s"^). 



3 PHYSICAL PARAMETERS AND 
DIAGNOSTIC DIAGRAMS 

Figure |2| shows the spatial profiles of low- to high-excitation 
oxygen emission lines ([O i] 6300A, [O 11] (3726A+3729A) 
and [O III] 5OO7A); we clearly see that the knots have a much 
lower excitation than the core. While the low-ionization pro- 
files ([O i] 6300A and [O 11] (3726A-H3729A)) show local 
maxima at the positions of the knots, the [O ill] 5007A pro- 
file presents its maxima in the core and does not peak at the 
knots. The average values of [N 11] 6583A/Hq at the knots 
are 2.0 (Knotl) and 2.1 (Knot2), much higher than those 
usually found for spherical and elliptical PNe (e.g., AUer & 
Czyzak 1983), although such values are frequently found in 
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the Peimbert Type I objects (Peimbert 1978; Peimbert & 
Torres-Peimbert 1983). The Core, on the other hand, has 
a less extreme [N ll] 6583A/HQf hne ratio, 1.2, even though 
still higher than in most PNe. Peculiar line ratios like those 
found in K 4-47 have been interpreted as indicative of the 
presence of shock-excitation or anomalous abundances (see, 
for instance, Miranda & Solf 1992). 

Fluxes were extinction-corrected using ch^i (the loga- 
rithmic ratio between observed and dereddened H/3 fluxes), 
determined from the observed Hq/H/3 and HQ/H7 ra- 
tios, for each of the selected features of the nebula, and 
the reddening law of Cardelli, Clayton, & Mathis (1989). 
The derived ch/3 are: 1.26±0.23(Knotl); 1.37±0.17(Core); 
l.ll±0.17(Knot2). The physical properties (for all the ne- 
bular regions) and chemical abundances of the Core were 
then evaluated assuming pure photoionization by the cen- 
tral star. 



3.1 Density and Temperatures 

We estimate the electron densities of K 4-47, which are listed 
in Table 1, using the [S 11] 6717A/6731A line ratio. Both 
knots are denser than the core (1900 cm"^) by factors of 2.4 
and 1.2 for Knotl and Knot2, respectively. 

Adopting these densities, electron temperatures are es- 
timated using standard line ratios: Te[0 ill] is obtained from 
the [O III] 4959A/4363A line ratio, which is appropriate to 
zones of medium to high excitation, and Te[N 11] from [N 
11] 6583A/5755A, representative of low-excitation regions. 
The Te[0 III] for Knotl and the re[N 11] of the Core, are 
lower limits. We find that, within the errors and when both 
determinations are available, the two temperatures are simi- 
lar. Errors on Tjfl ill] at the position of the knots are signif- 
icantly higher than those of the [N II] temperatures simply 
because of the low-ionization nature of these regions. Note 
that Te is remarkably higher than the typical values for PNe 
of around 10* K (Kaler 1986). 



3.2 Diagnostic Diagrams 

Figure 13 shows several diagnostic diagrams which are com- 
monly used to distinguish photoionized nebulae from shock- 
excited objects, adapted from Sabbadin, Minello & Bian- 
chini (1977) and Phillips & Cuesta (1999). In all of them, 
the knots of K 4-47 are located close to the shock-excited 
zones while the core is displaced toward the locus where 
PNe are usually found, i.e. photoionized objects. This 
suggests that shock-excitation plays an important role in 
the knots, but not in the core. The fact that in the [O 
III] (4959A+5007A)/H a vs. [S 11] (6716A-|-6731A)/Hq di- 
agram Knot2 do not fall exactly within the "Shock Regime" 
area is not against this interpretation, as its boundaries were 
originally computed for young stellar objects, using para- 
meters (gas abundances and shock velocities) which can be 
different from those of evolved objects like K 4-47. 

The large expansion velocities of the knots and their 
broadened line profiles (Corradi et al. 2000) are also typical 
of shocked fiows like H-H objects. In addition, the shock- 
excited nature of the knots would be consistent with their 
measured Te (Table 1), unusually large for low-ionization 
microstructures in PNe, and with their extremely high [N 
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Figure 3. Diagnostic diagrams showing the loci of the core and 
knots in K 4-47. They include the empirical loci of the PNe, super- 
nova remnants (SNR) and Hll regions (HIT), as well as the theo- 
retical "Shock Regime" zone as defined by Phillips and Cuesta 
(1999) using the plane-parallel and bow-shock models of Harti- 
gan, Raymond & Hartmann (1987). Line ratios considered are 
the sum of the doublets ([O lll] = [0 III] 4959A-I-5007A; [S ll] = [S 
11] 6716A-I-6731A; [N ll] = [N 11] 6548A-I-6583A) when not followed 
by the specific wavelength. 



11] 6583A/Ha line ratio. On the contrary, the velocities ob- 
served in the core are substantially lower than in the knots, 
and so is its [N 11] 6583A/Hq line ratio. 

It is therefore suggested that the high velocity knots 
Knotl and Knot2 of K 4-47 are mainly shock heated, while 
the Core is mainly photoionized. With the aid of existing 
photoionization and shock modelling we will explore in the 
following these alternatives. 



4 ANALYSIS OF THE CORE 
4.1 Empirical Abundances 

Ionic and total abundances for the Core were computed us- 
ing the ionization correction factors {icfj scheme, following 
Kingsburgh & Barlow (1994), as described in Corradi et al. 
(1997) and Goncjalves et al. (2003). Results are shown in 
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Table 2. Note that as the re[N ll] in Table 1 is only a lower 
limit, the abundances of the low-ionization ions are upper 
limits. Errors on the total abundances are obtained by pro- 
pagating the errors in the ionic abundances as well as on 
the icf, and are indicated within brackets in the table as 
percentage errors (14% for He, for which no icf is applied, 
and 30-50% for O and N). 

According to its large He and N/O abundances, K 4-47 
is similar to most bipolar PNe (Perinotto & Corradi 1998), 
and would be classified as a Type I FN (Pcimbcrt & Torrcs- 
Peimber 1983, see also Kingsburgh & Barlow 1994). How- 
ever, the O, Ne and S abundances of K 4-47 are among the 
lowest measured for Galactic PNc, and arc more typical of 
some halo (or Type IV) PNc (sec e.g. Howard, Henry & 
McCartney 1997). 

High He/H abundances (0.15 < He/H < 0.20), and some 
(but much less extreme than in K 4-47) oxygen deficiency 
(1.35 < O/H (xlO"") s; 2.7), were found in other bipo- 
lar PNe, such as NGC 6302 (Pottasch & Bcintcma 1999), 
NGC 6537 and He 2-111 (Pottasch, Beintema & Fcibclman 
2000). Very much like in the case of K 4-47, these PNe 
have considerably high electron temperatures (average Te, 
from optical, ISO and lUE spectra for the throe of them arc 
17200 K up to 18300 K). Therefore, we have the somewhat 
contradictory result that while K 4-47 seems to be an ex- 
treme bipolar PN in terms of its morphology and He and 
N abundances, its O, Ne and S abundances rather resemble 
PNe in the Galactic halo, as its relatively large height on the 
Galactic plane might also suggest (see Section 1), and con- 
trary to bipolar PNe which axe highly concentrated toward 
the Galactic plane (Corradi & Schwarz 1995). 

To further investigate these issues, in the next two Sec- 
tions wo will model the Core spectrum using, first, a shock 
excitation code, and, then, a pure photoionization model. 

4.2 Plane-parallel Shock Modelling 

Shock models have successfully reproduced the emission line 
ratios, line profiles and, recently, velocity maps in HH ob- 
jects (Beck et al. 2004). To check the possible influence of 
shock excitation on the Core spectrum we have explored a 
variety of steady plane-parallel models, for which we use the 
photoionization-shock code MAPPINGS Ic (Dopita, Binette 
& Tuohy 1984; Binette, Dopita & Tuohy 1985). 

One important characteristic of shock models is that 
the predicted spectra strongly depends on the preshock io- 
nization. We have compiled a set of fully pre-ionized shock 
models (i.e. the incident gas has been ionized) and a set of 
local equlibrium prcionization for the gas entering the shock 
(as described by Shull & McKee 1979, Hartigan et al. 1987). 
Relative abundances typical of Type I PNe were assumed, 
and two values for the pre-shock density were considered 
(100, 1000 cm~^). The input shock velocities were varied 
from 95 to 140 km s~^, i.e. the range producing large [O 
III] (4959A,5007A)/H/3 ratios (Hartigan et al. 1987), as ob- 
served. Note, however, that the measured expansion velocity 
of the Core is smaller than 50kms~'^ (Corradi et al. 2000), 
pointing to lower shock velocities. 

None of the predicted spectra could reproduce the large 
observed ratios of intermediate to high excitation emission 
lines with respect a H I Balmer line (i.e. [O II] 3727A/H/3, 
[O III] 5007A/H/3, [Ne iii] 3868A/H/3 and He ii 4686A/H/3). 



Moreover, the shocked spectra predicted by our models im- 
ply [O l] 6300 A/Ha and [S ll] (6716A-h673lA)/Ha larger 
than the observed values. A large discrepancy was also found 
with the observed [O ni] 4959A/4363A emission line ratio, 
which is more than 2 times larger than the model values. 
And, finally, [N ii] 6583A/5755A is very high, « 300, be- 
cause the intensity of [N ii] 5755A is largely understimated 
by these models. 

We conclude that the plane-parallel shock models axe 
not able to reproduce the spectrum of the core of K 4-47. 

4.3 Photoionization Modelling 

We have used the photoionization code CLOUDY 95.06 
(Ferland et al. 1998). 

CLOUDY needs as input the information on the shape 
and intensity of the radiation from the ionizing source, the 
chemical composition and geometry of the nebula, as well as 
its density and size. As mentioned in Section 1, the distance 
of K 4-47, and thus its size and the luminosity of the central 
star, are poorly known. We adopt the distance of 5.9 kpc 
computed by Tajitsu & Tamura (1998), although extending 
the calculations to the full range of distances from 3 kpc to 
7 kpc proposed by Corradi et al. (2000). 

A spherical geometry, and a filling factor of 1.0 have 
been assumed. Using the Hq and [N ll] images from Corradi 
et al. (2000) a core diameter of 1.9 arcsec was estimated 
(10% contour extension corrected for the finite resolution; 
Tylenda et al. 2003). The density was initially kept to con- 
stant and equal to its empirical value (1900 cm~^, see Ta- 
ble 1). However, as noted in the Introduction, Aaquist & 
Kwok (1990) discovered a bright and very compact radio 
core at the centre of K 4-47, implying a high density of 
72000 cm~^ (assuming an optically thin nebula at 5 GHz, 
the distance of 5.9 kpc and following Goudis 1982). For this 
reason, models with much higher density values were also 
explored. 

Dust grains have been included, since they have an im- 
portant effect, particularly on the temperature structure of 
PNe. This effect depends on the type of grain (graphite 
and/or silicate), the grain abundances, and the grain-size 
distribution, being more relevant in the inner regions of the 
nebula (see Dopita & Sutherland 2000). Both graphite and 
silicate grains were considered. As we do not know the grain 
size distribution and grain abundances for K 4-47 and more 
in general for PNe, we adotped ISM size distribution, with 
ISM gas-phase depletions due to grains, following van Hoof 
et al. (2004). 

We adopted a blackbody spectrum for the central 
star, with Te// derived from the H/3, Hel and Hell emis- 
sion line ratios, using a modified Zanstra method as done 
for instance by Mikolajewska et al. (1999). This gives 
Te//=130 000 K (from Hell 4686A/H/3) and re//=115 000 K 
(from Hel 4686A/5876A). A lower hmit to the central star 
luminosity comes from the IRAS spectral energy distribu- 
tion (Tajitsu & Tamura 1998), which gives 16 D^ L© 
(where D is the distance in kpc) yielding 550 Lq at 5.9 kpc. 

As with the elemental abundances, we first run models 
with the empirical abundances given in Table 2, and sub- 
sequently with "average" values for either normal (Type- 
II), or Type-I PNe, following Kingsburgh & Barlow (1994). 
CLOUDY models with the empirical abundances do not 



K J^-J^l: photons versus shocks 5 



match the observed hne fluxes. Both hehum and [N 
ii] (6548A,6583A) emission lines are largely overestimated, 
while oxygen, sulphur and neon lines are underestimated, 
mainly because the very low 0/H, Ne/H and S/H input 
abundances arc not compensated by a high electron tem- 
perature in the model as empirically determined. Tests with 
Type-II PNe abundances were also not successful because 
the predicted [N ii] 6548A,6583A intensities are very low, 
which is expected as the main difference between Type-I 
and Type-II abundances is that the latter arc depleted in 
nitrogen by a factor of ^4 (Kingsburgh & Barlow 1994). 
Assuming Type-I abundances, we found a partial agree- 
ment with the observed spectrum for r,.// = 120 000 K, 
L=550 Lq and Ne= 1900 cm^''. Important optical lines like 
Hen4686A, [O iii] (4959A.5007A), [N ii] (6548A,6584A) and 
[S ii] (6716A,6731A) are well reproduced by this model, with 
discrepancies of 10% or less. Hei 447lA, 6678A, Hell 5412A, 
and [S III] 6312A lines agree with the observed spectrum 
within 12% to 20%. Other relatively intense lines such as 
[O II] (3726A+3729A), [Neiii] 3868A, and Hei 5876A show 
discrepancies of 35% to 50% with respect to the model. Fur- 
thermore, the model underestimates the intensities of the [O 
III] 4363A and the [N ll] 5755A lines by a factor of about 3. 
These lines arc the keys for the determination of the electron 
temperature. Finally, the model underestimates [Nl] 5200A 
and [Oi] 6300A by factors larger than 30. Let us note, in ad- 
dition, that this model requires a size for the core of 4.0 arc- 
sec, twice the optical size observed. 

It is then instructive to explore alternative scenarios. 
Let us first note that in some PNe and related objects (see 
e.g. Corradi 1995), high [O ill] 4363A and [N ll] 5755A line 
intensities have been interpreted as the signature of very 
high core densities, as at A^e larger than about 10'' cm^'^ 
the auroral ([O iii] 4363A and [N ii] 5755 A) to nebular ([O 
in] 4959A and [N 11] 6583A) line ratios are indicators of 
density rather than of temperature (Gurzadyan 1970). For 
this reason, we have calculated other models assuming much 
higher (albeit still constant) densities. As noted above, high 
core densities are implied by the radio size and flux (Aaquist 
& Kwok 1990). 

CLOUDY models for the observed radio core size, 
0.25 arcsec, and such high densities (from 72000 K up to 
3.0x lO^cm"^) show that both [O ill] 4363A and [N 11] 5755A 
intensities can be now reproduced, but other nebular lines, in 
particular [O 11] (3726A-h3729A), [N 11] (6548A,6583A), and 
[S 11] (6717A,673lA), become now largely underestimated 
because of coUisional quenching. A natural way to solve the 
problem might be to assume a strong density stratification 
in the core of K 4-47, with a very dense inner zone where 
[O III] 4363A and [N 11] 5755A are mostly formed, and a 
lower density outer region where other important nebular 
lines axe produced. This idea should be tested by means of 
an appropriate photoionization model; 3-D codes (like MO- 
CASSIN; Ercolano et al. 2003) are much better suited than 
CLOUDY to deal with such extreme density variations. 

In summary, we find that none of the constant density 
models is able to account, simultanoouly, for all optical emis- 
sion lines in the Core. In particular, the [O ill] 4363A and 
[N II] 5755A intensities arc strongly understimatcd if the 
nebular density is the one derived empirically from the [S 
11] lines. A model with a strong density stratification could 
possibly offer a solution to the problem. 
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Figure 4. Results of the bow-shock models, in terms of line ra- 
tios, as compared to the observed emission. Results are shown 
for bow-shock velocities of Ybs of 50, 75, 100, 150, 200, 250, and 
300 km s^ ^ , corresponding to points moving from left to the right 
side inside each panel. 

5 SHOCK MODELLING FOR THE KNOTS 

The fact that we are not able to provide an accurate repre- 
sentation of the Core spectrum, also prevents us to deter- 
mine the radiation escaping from the core and reaching the 
knots, and thus to attempt a reliable photoionization mo- 
delling of the latter. A consistent photoionization modelling 
of the whole nebula (core-|-knots) would also require a fully 
3-D modelling, that goes beyond the scopes of the present 
work. Given the evidence that the knots might be excited by 
shocks (Section. 3.2), we instead attempt to describe their 
spectrum using existing shock models. 

5.1 Bow-Shock Models 

We have used the bow-shock models described in Raga & 
Bohm (1986) and Hartigan et al. (1987). We consider a bow 
shock with a functional form, z/a = {r/aY, where z is mea- 
sured along the symmetry axis, r is the cylindrical radius, 
the parameter p determines the shape of the bow shock and 
the constant a determines its size (Beck et al. 2004). The 
emission from the bow-shock is modeled as described by 
Hartigan et al. (1987), being the geometry of the bow shock 
the parameter that determines the shock velocity. 

The plane-paralel shock models required to predict 
the bow-shock emission line ratios are obtained with 
the photoionization-shock code MAPPINGS Ic. We have 
adopted a prc-shock density of 400 cm~^, that is within 
the range generally assumed for stellar jets. The pre-shock 
magnetic field was taken to be negligible. We assumed ele- 
mental abundances for He, C, N, O, Ne, S, and Ar to be 
within the range of Type-I PN values. As for the preionizar 
tion of the gas, local equilibrium is assumed. We consider 
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the bow-shock velocity as a free parameter, adopting a set 
of values from 50 to 300 km s~^. 

5.2 Comparison with the Observed Spectra 

We have obtained the best fit to the observed spectra for 
p=3. A summary of the results is presented in Figure |1| for 
a number of emission line ratios plotted as a function of 
[N ll] 6583A/Ha. Note that the later line ratio increases in 
strength as the shock velocity increases, and in the case of 
K 4-47 is matched by only the largest shock velocities con- 
sidered, somewhat larger than the estimate of 150 kms~^ 
by Corradi et al. (2000). 

For these velocities, the bow-shock models fit pretty well 
the line ratios [S ii] (4069A-H4076A)/(6717A+673lA)=[S 

II] 4072A/6725A, and [N ll] 6583A/5755A, sensitive to the 
temperature in the recombination zone where [S ll] and [N ll] 
are formed, but not equally well the [O ill] 4959A/4363A ra- 
tio, indicating that the predicted electron temperature re[0 

III] is too large. The [N ll] and [O ill] electron temperatures 
predicted by these models, for both knots, are 17800 K and 
27600 K, respectively. 

The ratios [O l] 6300A/Hq and [O 

II] (3726A+3729A)/H/3 are also well reproduced. Larger 
discrepancies appear for the helium lines Hel 5876A and 
Hen 4686A, and specially for [Ne ill] 3868A and [O 

III] 5007A in Knot2, and [S ll] (6716A-^673lA) and [N 
l] 5200A in both knots. Note that a similar problem for 
the [S ll] (6717A-h673lA)/Ha ratio is also found in HH 
objects (Raga et al. 1996). The large discrepancy for the [N 
l] 5200A/H/3 emission line ratio is intriguing, as strong line 
intensities for this ion can only be produced with low shock 
velocities, which could however not account for the observed 
fiux of higher ionization stages like [N ll] 6583A. Diff'erent 
shock geometries, producing a different distribution of 
shock velocities, might help to solve this problem. 

Also note that the bow-shock models are based on ap- 
praisal of what would be reasonable parameters for the knots 
of K 4-47, and we did not intend to reproduce the line inten- 
sities separately for each of the knots. But the knots are not 
identical: in particular, high excitation lines like Hell 4686A, 
[Ne III] 3868A and [O lll] 5007A are much brighter in Knot2 
than in Knotl, which indicates that the radiation field of 
the central star would affect the Knot2 emission line spec- 
trum much more than that of Knotl. The observed spectra 
of Knot2 would also fit a model with a somewhat larger N/H 
abundance and a lower bow shock velocity. 

Summarising, the analysis of the shock-excited knots 
suggests that they are high-velocity (250 to 300kms~^) 
condensations with high electron temperatures (TefN ll] 
= 17800 K) and chemical abundances similar to Type-I PNe. 



6 DISCUSSION 

The analysis presented in this paper is the first attempt 
to determine the physical parameters, excitation mecha- 
nism and chemical abundances of K 4-47. Among other PNe 
with LISs, M 2-48 (Lopez-Martin et al. 2002) and Kj Pn 8 
(Lopez, Vazquez & Rodriguez 1995) have, similarly to K 4- 
47, a high electron temperature and the evidence for shock- 
excited LISs. But note that in these PNe, except for one of 



the M 2-48's knots (for which Te is 20100 K, twice the core 
temperature), Te and Ne are not estimated in a detailed 
zone-by-zone basis, because [O ill] 4363A and [N ll] 5755A 
were not measured. 

The two distinct models that we used for the analysis 
of K 4-47 (pure photoionization for the core and pure shock- 
excitation for the pair of knots), provide some clues for the 
interpretation of this nebula. 

First, the analysis of the core would suggests He and N 
overabundances typical of Type-I (Galactic disk) PNe, but 
also an extreme oxygen deficiency which is instead more ty- 
pical of PNe in the Galactic halo. These abundances calcula- 
tion would however be invalidated if more detailed modelling 
proves the existence of a strong density stratification within 
the core, with an inner region with extremely high density 
(~ 10^ cm~^), and an outer zone with a much lower density 
matching the empirical determination from the [S ll] lines. In 
fact, in this case the Te determined empirically using the [O 
III] 4959A/4363A and [N ll] 6583 A/5755 A line ratios would 
be wrong, as the auroral and nebular lines involved would be 
formed in different regions. Note that, albeit unusual, such 
high density cores were actually found in several other PNe 
(Corradi 1995), which like K 4-47 show an extreme degree 
of coUimation. 

Second, the knots seem to be mainly shock-excited. This 
conclusion is strenghtened by diagnostic line ratio diagrams, 
and, interestingly, by the fact that the H2 emission of K 4-47 
is also shock-excited, contrary to most PNe where instead it 
is excited by fluorescence (Lumsden et al. 2001). Our shock 
modelling also indicates that knots have Type-I abundances. 
According to these models, knots would move with space 
velocities of 250-300 kms~^, about twice the value estimated 
by Corradi et al. (2000) from the analysis of the lines profiles. 
Moreover, while the shock model is successful in matching 
the [N 11] (6583A/5755A) line ratio at the knots it underes- 
timate their [O ill] (4959A/4363A) ratio by some 30% and 
70% for Knotl and Knot2, respectively. 

As a final remark, the two physical processes explored 
separately in this paper, i.e. photoionization and shock heat- 
ing, are likely to be simultaneously present in a real nebula 
like K 4-47, as gas can be ionized by both energetic pho- 
tons from the central star and from the shocks associated 
with the observed supersonic outflows. Taking this into ac- 
count would require a sophisticated modelling which is not 
presently available, and which would also need to be sup- 
ported by a better knowledge of the nebular and stellar pa- 
rameters. In particular, very little information is presently 
known about the distance and the properties of the central 
star of K 4-47, a basic data that could be addressed by future 
observations. 
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Table 1. Observed emission line fluxes (normalized to H/3 
=100), electron density and temperatures. The estimated per- 
centage errors in the line fluxes (for a range of fluxes relative 
to H/3) are given at the bottom. 



Line Fluxes 



Line Identification (A) 


Knotl 


Core 


Knot2 


[On] 3726.0 + 3728.8 


23.4 


26.6 


33.3 


H9 3835.4 


9.46 


- 


- 


[Nciii] 3868.7 




48.9 


40.3 


Hcl 3888.7 + H8 3889.1 




7.05 


12.8 


[Nciii] 3967.5 + He 3970.1 


11.5 


38.9 


27.8 


[Sll] 4068.6 


23.1 


3.18 


20.0 


[Sll] 4076.4 


15.0 


2.49 


10.7 


US 4101.8 


15.6 


16.1 


16.6 


H7 4340.5 


30.1 


32.0 


27.8 


[Oiii] 4363.2 


3.45 


20.8 


4.56 


Hel 4471.5 




4.28 


9.22 


Niii + On 4641.0 


3.69 


4.40 




Hen 4685.7 


2.95 


24.4 


14.3 


H/3 4861.3 


100 


100 


100 


[Oiii] 4958.9 


35.9 


361 


100 


[Oiii] 5006.86 


118 


1088 


279 


[Ni] 5200.2 


367 


42.0 


421 


Hell 5411.6 




4.11 


- 


[Oi] 5577.4 


12.6 




10.4 


[Nil] 5754.6 


51.4 


43.6 


40.3 


Hcl 5875.7 


11.5 


42.5 


14.6 


Hcll 5882.0 






5.54 


Hen 5896.5 


5.82 




10.1 


[Oi] 6300.3 


427 


54.7 


305 


[Sm] 6312.1 




7.65 


- 


[Oi] 6363.8 


146 


17.1 


99.1 


[Nil] 6548.0 


446 


333 


415.2 




681 


762 




[Nil] 6583.4 


1406 


976 


1276 


Hoi 6678.1 


3.22 


13.4 




[Sll] 6716.5 


93.5 


32.1 


86.0 


[Sll] 6730.8 


165 


45.2 


130 


Ne[Sll](cm-3) 


4600±850 


1900±420 


2400±400 


Te[OlIl](K) 


^21000 


19300±2300 


16100±4400 


Te [Nll](K) 


18900±2950 


^21000 


16950±2800 




Percentage errors in line 


fluxes 


Line Fluxes 


Knotl 


Core 


Knot2 


(0.01-0.06)Ih3 


51 


41 


58 


(0.06-0.15)Ih3 


23 


15 


26 


(0.15-0.30)Ih/3 


14 


11 


15 


(().30-2.0)Ih3 


12 


8 


11 


(2.0-5.())1h3 


12 


7 


9 


(5.0-10.0)Ih3 


11 


7 


9 


> 10 Ih/3 


11 


7 


9 
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Table 2. Ionic and total abun- 
dances of the Core. Percentage 

errors arc gi\'on within brarkcts. 



Elemen /lont 


Abundances 


Hc+ /H 


1.14E-01(08) 


Hc2 + /H 


2.56E-02(14) 


He/H 


1.39E-01(14) 


O^/H 


6.10E-06(10) 


0+/H 


3.09E-06(14) 


0^+/H 


6.13E-05(21) 


ic/(0) 


1.14 


O/H 


r.37E-05(32) 


N^/H 


8.47E-06(14) 


N+/H 


l.B7E-05(8) 


jc/(N) 


23.8 


N/H 


3.74E-04(40) 


Ne^+/H 


1.44E-05(32) 


ic/(Ne) 


1.20 


Ne/H 


1.74E-05(66) 


S+/H 


2.24E-07(14) 


S^+/H 


7.44E-07(36) 


ic/(S) 


2.02 


S/H 


1.96E-06(48) 



